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Abstract. Pelophylax esculentus is a vertebrate animal with hemiclonal heredity, attracting the interest of many reproduc¬ 
tive and evolutionary biologists. It is a hybrid between P. ridibundus and P. lessonae. These three taxa form the so-called 
Pelophylax esculentus complex with a population structure usually comprising a hybrid taxa and one parental species. Data 
on population types at the southernmost distribution area of their sympatry are rare. Fiere we sampled five sites in inland 
Croatia in order to analyse the population structure, sex ratio and age structure. The individual genotypes of 93 randomly 
collected water frogs were verified with allozyme markers for three species-specific polymorphic loci. In order to estimate 
population age structure, the annual growth rate (skeletochronology) and growth index profiles were also investigated. The 
growth index profiles were analysed by an estimation of number of lines of arrested growth visible in the cross-section of 
femur bones. Our results revealed the presence of the R-E-L population with a dominance of P. esculentus. Pelophylax ridi¬ 
bundus was the least abundant taxon but with a relatively high age estimate of eight years on average. Its annual growth 
rate did not differ from the remaining two species. Gene introgression of mostly ridibundus alleles was also observed in 
hybrids. Most profiles of gonads in hybrids showed presence of both parental genomes with dominance of ridibundus al¬ 
leles. The study area represents one of the southernmost distributions of the hybrid taxon in Europe, making it attractive 
to study gene flow and impact of P. esculentus on P. ridibundus, a typical water frog representative of the Balkan Peninsula. 

Key words. Hybridogenesis, population composition, allozyme markers, hemiclonal reproduction, gene introgression, 
skeletochronology. 


Introduction 

The Edible Frog, Pelophylax esculentus (Linnaeus, 1758), is 
of particular interest to evolutionary biologists and popu¬ 
lation ecologists because of its hybrid origin and unusual 
reproduction mode. It represents a natural hybrid taxon 
between the Pool Frog, P. lessonae (Camerano, 1882), and 
the Marsh Frog, P. ridibundus (Pallas, 1771). The genom¬ 
ic constitutions of the parental species are denoted as LL 
and RR, respectively. Pelophylax esculentus reproduces via 
hybridogenesis. In this hemiclonal reproduction, a clonal 
copy of either lessonae or ridibundus genome is transmit¬ 
ted to gametes (Graf & Polls-Pelaz 1989, Bergen et al. 
1997, Holenweg Peter et al. 2002, Pruvost et al. 2013). 
The other genome is excluded during the first division of 


gametogenesis (Tunner & Heppich-Tunner 1991). There¬ 
fore, hybrids have to backcross with a parental species to 
receive an excluded genome and produce a new generation 
of hybrids (Graf 1986, Semlitsch et al. 1996, Vorburger 
& Reyer 2003, Ghristiansen & Reyer 2011). As a result, 
P. esculentus maintains a permanent Fi genotypic constitu¬ 
tion. 

The compositions of population systems in water frogs 
vary throughout Europe (Plotner et al. 1994, Holsbeek 
& Jooris 2009, Ghristiansen & Reyer 2011, Mayer et al. 
2013, Pruvost et al. 2013, Hoffmann et al. 2015). The type 
of population system depends on the environmental var¬ 
iables and genetic background of the hybrids (Sas 2010). 
Pelophylax esculentus usually has a diploid constitution and 
persists in mixed populations with P. lessonae (here called 
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the L-E system) or P. ridibundus (the R-E system) (Graf & 
Polls-Pelaz 1989). Alternative population systems such as 
those that involve the presence of triploid individuals (LLR 
and ERR) forming all-hybrid populations (E-E system) are 
also known, mostly from northwestern Europe (Chris¬ 
tiansen & Reyer 2011, Pruvost et al. 2013, Pruvost et al. 
2015, Hoffmann et al. 2015). The L-E system represents the 
most widespread population type (Graf & Polls-Pelaz 
1989). Relatively rare are populations where all three taxa 
coexist within the so-called R-E-L system found in Latvia 
(Borkin et al. 1986), Switzerland (Hotz et al. 1992), Hun¬ 
gary (Gubanyi 1992), central and northwestern Ukraine 
(Mezhzherin & Morozov-Leonov 1993), Czech Re¬ 
public (Kotlik & SuLOVA 1994), European Russia (Lada 
et al. 1995, Borkin et al. 1986, 2002), western Germany 
(ScHROER & Greven 1998), Romania (Sas 2010), Slovakia 
(Pruvost et al. 2013, Mikulicek et al. 2014), northern Ser¬ 
bia (Spasic-Boskovic et al. 1999, Krizmanic & Ivanovic 
2010) and Croatia (Berger et al. 1988). Pelophylax ridibun¬ 
dus is assumed to be ethologically isolated in the R-E-L 
system (Hotz et al. 1992). It is known that in some cases 
P. ridibundus migrates toward ponds with P. lessonae and 
P. esculentus during the breeding period converting the L-E 
system into the R-E-L system (Sas 2010). 

After the initial recognition that P. esculentus is of a hy¬ 
brid origin (Berger 1966,1967), a number of studies were 
performed on morphological investigations, experimental 
hybridizations (e.g. Berger 1968, Gunther 1973, Eben- 
DAL 1979), cytological analyses (e.g. Heppich & Tuner 
1979, Heppich et al. 1982) and molecular genetic analy¬ 
ses (e.g. Engelman 1972, Graf & Muller 1979, Graf 
& Polls-Pelaz 1989, Hauswaldt et al. 2012) mostly on 
populations from central, western and northwestern Eu¬ 
rope (Plotner 2005, Pruvost et al. 2015 and citations 
therein). 

Considering P. esculentus hybrid's ability to serve as a 
vehicle for the transfer of P. lessonae genes into P. ridibun¬ 
dus in the area of sympatry (Plotner et al. 2008), suggest¬ 
ing the impact of hybrids on the parental taxon especially 
north from the Balkan Peninsula, our knowledge on P. es¬ 
culentus distribution is important. Previous studies sug¬ 
gested that water frogs are indigenous in Croatia (Jelic et 
al. 2012, Jelic 2014) and morphological analyses point out 
a presence of all three taxa in the region (Gasc et al. 1997, 
Jelic et al. 2012, Karaica et al. 2016). In this investigation 
we aimed to describe the population structure of the P. es¬ 
culentus complex in the inland part of Croatia from dif¬ 
ferent aspects. Its taxonomic composition was investigated 
by means of allozyme analyses, whereas skeletochronology 
was employed in order to investigate the age structure of 
the complex population in question. Furthermore, we have 
analysed the sex structure of this population. All investiga¬ 
tion aspects: genetic, age and sex structure, are important 
in describing a status of a certain population, in particular 
its ecological relationships and stability. Since population 
structure and demographic parameters have not been ana¬ 
lysed for populations in Croatia, the unexplored southern 
limits of this complex’ range, we believe that obtained data 


could contribute significantly not only in widening knowl¬ 
edge on this unusual biodiversity component, but also for 
making scientific bases for its conservation. 

Materials and methods 

Sampling 

The sample sites have been described in detail by Karai¬ 
ca et al. (2016). In brief, the sampling took place near 
the Ilova River (Danube River system) in the northwest¬ 
ern part of Croatia (Bjelovar-Bilogora County; 100-125 m 
a.s.l.) at five localities described as follows: the pond near 
the fish hatchery (Crubisno polje: 45°42’3i” N i7°o8’52”E), 
two sites located near the Ilova River channel (Ulovcev 
mlin: 45 ° 39 ’ 54 ”N, 17 ° 13 ’ 36 ”E; Mali Zdenci: 45°4o’io”N, 
17°o8’26”E), and the wet meadows near the brook (Ve- 
lika Barna: 45°44’o9”N, i7°o6’53”E; Velika Jasenovaca: 
45°43’4i”N, 17°o6’oi”E). Distances between the localities 
ranged from 4.5 to 8 km. The collection of samples was 
approved and carried out with the permission of the Di¬ 
rectorate for Nature Protection, Ministry of Culture of the 
Republic of Croatia (permit number given in the acknowl¬ 
edgements). Frogs were collected by a hand net during day 
or night. Sex was determined by the presence or absence of 
vocal sac openings, thumb pads and/or by direct inspec¬ 
tion of gonads. Those frogs that had not yet developed the 
external sex characters and, simultaneously, their inspec¬ 
tion of gonads were not done, we further call as subadults. 
Taxon identification was first based on external morpho¬ 
logical characters (Gunther 1990, Plotner 2005). The re¬ 
sults of our morphology investigation were in details de¬ 
scribed in Karaica et al. (2016). The determination of the 
species was later verified by isozyme (allozyme) electro¬ 
phoresis. Samples of somatic and gonadal tissues were iso¬ 
lated and frozen at -20°C for further genetic observation. 
A Chi-square test was used in order to detect possible dif¬ 
ferences between the observed population sex ratio and the 
expected 1:1 sex ratio. Software Packages Microsoft® Excel 
2010 and StatSoft® Statistica Version 7 were employed for 
data analyses. The frogs were marked and preserved in 75% 
ethanol at the Department of Zoology, Faculty of Science, 
University of Zagreb. 

Allozyme analysis 

Approximately icm^ of skeletal muscles and gonads were 
crushed separately and homogenised in an equal volume 
of pH 8.5 Tris NaCl extraction buffer (Valenta et al. 1971, 
Appendix II). Following cold centrifugation at 4°C, equal 
volumes of the supernatant were mixed. The mixture was 
frozen overnight in a deep freezer at -70°C. The enzymes 
were separated horizontally using gel electrophoresis 
method in a refrigerator at 4°C. The enzymes were separat¬ 
ed on 11% starch gel using a pH 8.6 Tris-citrate buffer sys¬ 
tem (Valenta et al. 1971, Appendix II) and pH 6.0 Tris-ci- 
trate acid electrode buffer (Uzzell & Berger 1975, Appen- 
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dix II). The supernatants mixtures were applied to the gels 
on filter paper tabs (Whatman no 3.). The investigated loci 
were already described as polymorphic and species-specif¬ 
ic (UzzELL & Berger 1975, Beerli 1994) for the following 
enzymes: aspartate aminotransferase (Aat; EC 2.6.1.1), lac¬ 
tate dehydrogenase (Tdh-i; EC 1.1.1.27) and phosphoglu- 
comutase (Pgm-2; EC 5.4.2.2). The enzymes were stained 
using procedures similar to those described by Harris & 
Hopkinson (1976), Buth & Murphy (1980) and Pasteur 
et al. (1987, Appendix II). We used a previously determined 
water frog as a biological standard for allele mobility at 
every allozyme run. The allele homology was identified af¬ 
ter KotlIk & SuLOVA (1994) for Aat, Uzzell et al. (1980) 
for Pgm-2, and Vorburger (2001) and Vorburger & Ul- 
rich-Reyer (2003) for Ldh-i. In our analysis, we assumed 
that bands with equivalent mobility are controlled by the 
same allele. Stained gels were photographed. The fastest 
visualised allele products were designated as a. Samples 
that revealed unclear patterns were reprocessed again. The 
relative intensity of bands has also been studied to estimate 
how many copies of the gene coding the product were ex¬ 
pressed in the tissue of each individual. Here, the asym¬ 
metrical staining intensity would potentially correspond to 
triploids (Tunner 2000). Allele frequency was measured 
as the proportion of a particular allele among all of the al¬ 
lele copies being considered. 

Skeletochronology 

In order to estimate the age of the 93 genotyped water frogs 
from our population sample, we used a skeletochronolo¬ 
gy method that has already been applied on a variety of 
amphibian species from temperate and tropical regions 
(Kumbar & Pancharatna 2001, Lin & Hou 2002, Yil- 
MAZ et al. 2005, Guarino et al. 2008). During dissection, 
we successfully isolated 81 femoral bones without any dam¬ 
age, while the remaining 12 bones had damaged diaphy- 
sis and were excluded from the analysis. The bones were 
stored in 70% ethanol for further processing. Each bone 
was then transferred into distilled water for 24 h, fixed 
with 10% formalin for 24-48 h and decalcified with Rapid 
Decalcifier solution (Apex Engineering Product Corpora¬ 
tion) for 5-8 h. After decalciftcation, the residual decalci¬ 
fying solution was removed under running tap water. The 
processed bone tissue was cross-sectioned at a thickness 
of 20-25 pm at -25°C using a Leica Cryocut 1800 cryostat. 
These cross-sections were stained with Ehrlichs hematox¬ 
ylin in the dark for 5-10 min, rinsed in water two times 
for another 5-10 min, then mounted in a drop of glycerine 
gelatine on glass microscope slides and covered with cover 
slides. Sections were examined with a Zeiss AxioVert 200 
microscope (Carl Zeiss Microimaging) and photographed 
with a mounted Zeiss Axiocam MRC camera using Carl 
Zeiss Axiovision Rel.4.6 software (Carl Zeiss Microimag¬ 
ing GmbH). 

Estimations of individuals’ age were done on at least 
four mid-diaphysis sections which have the thickest cor¬ 


tical bone. The assumption was that the number of bright 
zones followed by so-called dark lines of arrested growth 
(LAGs; developed during hibernation period) seen in the 
periosteal bone represents growth during a single year 
(Castanet & Smirina 1990, Erismis 2011). Lines of arrest¬ 
ed growths on bone sections were difficult to read in two 
individuals (one from P. ridibundus and one from P. escu¬ 
lentus), thus they were not included in the statistical analy¬ 
sis. While examining the photographs of the bone cross- 
sections, two authors counted the LAGs independently. 
When discrepancies in the estimated number of LAGs oc¬ 
curred, the authors rechecked photographs of that partic¬ 
ular cross-section and came to a consensus regarding the 
final number. Nonparametric statistical Kruskal-Wallis test 
(a = 0.05) was used to detect possible differences in esti¬ 
mated age (age structure) and annual growth rate data (to¬ 
tal body length/estimated age ratio) between all three wa¬ 
ter frog species in the studied Pelophylax complex popula¬ 
tion. 


Results 

Population composition 

All three taxa, P. lessonae, P ridibundus and P. esculentus of 
both sexes, were found in the Ilova River drainage, Croatia. 
Out of the 109 specimens caught, 93 (85.3%) in total were 
electrophoretically determined with species-speciftc allo¬ 
zyme markers. The remaining 16 specimens (14.7%) were 
not analysed due to a lack of tissue for analysis or inactiv¬ 
ity of somatic tissue enzymes. Using the allozyme markers, 
the majority of specimens (66.7%) were identified as P. es¬ 
culentus, while 25.8% belonged to P lessonae and 7.5% to 
P. ridibundus (Fig. 1). 

Out of the 93 frogs, 61 were adult frogs from which 26 
were females (27.9% of totally collected frogs) and 35 were 
males (37.6% of totally collected frogs). There were 32 sub¬ 
adults where sex determination was impossible (34.4%). 
Female to male sex ratios did not differ significantly from 
the expected 1:1 ratio (Chi-square, df = 1, p = 0.25). Among 
the collected P. esculentus (n = 62) individuals, 27.4% were 
females (n = 17), 37.1% were males (n = 23) and 35.5% were 
subadults (n = 22). Out of the 24 P. lessonae individuals we 



P. esculentus P. lessonae P. ridibundus 


Figure 1. Sex composition and proportion of subadult individuals 
of the studied Pelophylax individuals in Croatia. Abbreviations: 
F - females, M - males. Sub - subadults. 
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Table 1. The allele frequencies investigated at the three loci for Pelophylax ridibundus, P. lessonae and P. esculentus. Abbreviations: 
Aat - aspartate aminotransferase, Ldh-1 - lactate dehydrogenase, Pgm-2 - phosphoglucomutase. 


Locus 

Aat 

Pgm-2 


Ldh-1 


Allele 

a 

b 

b c 

a 

b 

c 

d 

P. ridibundus 

1.000 

0.000 

0.000 1.000 

0.250 

0.000 

0.750 

0.000 

P. lessonae 

0.000 

1.000 

1.000 0.000 

0.000 

0.429 

0.000 

0.571 

P. esculentus 

0.510 

0.490 

0.578 0.422 

0.096 

0.061 

0.588 

0.254 


identified, 29.2% were females (n = 7), 37.5% males (n = 9), 
and 33.3% subadults (n = 8). Pelophylax ridibundus (n = 7) 
consisted of 28.6% females (n = 2), 42.8% males (n = 3), and 
28.6% subadult individuals (n = 2) (Fig. 1). 

Based on the pattern of the species-specific allozyme 
loci, three P. ridibundus males, six P. lessonae (four males, 
two females) and 29 P. esculentus (14 males, 9 females and 
six subadults) were sampled in Grubisno polje, 16 P. lesso¬ 
nae (five males, three females, and eight subadults) and 25 
P. esculentus (three males, seven females and 15 subadults) 
were sampled at Ulovcev mlin. Two subadult P. ridibun¬ 
dus, two P. lessonae females and seven P. esculentus (six 
males and one female) were found at the locality of Mali 
Zdenci. A single subadult P. esculentus was found at Velika 
Barna and two P. ridibundus females were found at Velika 
Jasenovaca. 

Allozyme genotype composition 

Aat. Two alleles typical for both species were found (a and 
b). All of the P. ridibundus were homozygous for the allele 
a, while all of the P. lessonae were homozygous for the al¬ 
lele b. Most of the P. esculentus individuals showed a het¬ 
erozygous combination (a/b), while one individual had a 
homozygous combination typical for P. ridibundus (a/a). 
The staining intensity of bands did not appear to differ be¬ 
tween alleles. 

Pgm-2. Two alleles (b and c) were found. Parental spe¬ 
cies showed homozygous allele combinations for this lo¬ 
cus. Pelophylax lessonae individuals expressed the faster 
moving allele b, while P. ridibundus individuals showed the 
expression of the slower allele c. Hybrids showed hetero¬ 
zygosity in most cases, but there were 9 hybrids express¬ 
ing alleles typical for P. lessonae (b/b) individuals, while 
only two contained a homozygous combination typical 
for P. ridibundus (c/c). The staining intensity of bands was 
symmetrical. 

Ldh-i. At this locus, four alleles were found (a, b, c and 
d). The investigated P. ridibundus individuals expressed c/c 
and a/c genotypes, while P. lessonae had b/b, d/d or b/d 
allele combinations. Hybrids showed all kinds of varia¬ 
tions of heterozygosity between the two parental species 
(a/d, b/c and c/d), except in two individuals that had an 
expression of a/c alleles and 14 individuals that had an ex¬ 
pression of c/c alleles typical for P. ridibundus. The staining 
intensity of bands was asymmetrical in 14 hybrid individu- 


Table 2. Summary of detected allele gene introgressions in Pelo¬ 
phylax esculentus. The underlined letters represent aUeles specific 
for P. ridibundus-, non-underlined aUeles typical for P. lessonae. 
Abbreviations: N - number of individuals; Aat - aspartate ami¬ 
notransferase, Pgm-2 - phosphoglucomutase, Ldh-1 - lactate 
dehydrogenase. 


N 

Aat 

Pgm-2 Ldh-1 

1 

a/a 


9 


b/b 

2 


c/c 

14 


c/c 

2 


a/c 

1 


a/c 


als, of which one frog showed b/c, one a/c, four a/d and 
eight individuals expressed the c/d allele combination pat¬ 
tern (more intensely stained bands on the gel electrophore¬ 
sis are in bold). This phenomenon was observed in 34.8% 
males and 11.8% females. 

The allele frequencies for each taxon are summarized in 
Table 1. 

The gene introgression (gene flow from one taxon to an¬ 
other) was found in 24 hybrid individuals among the three 
loci examined (38.7% of all hybrids). The mostly intro- 
gressed was the Ldh-i locus with 17 individuals expressing 
P. ridibundus-like alleles only. There were 14 hybrids ho¬ 
mozygous for the allele c (Grubisno polje, Ulovcev mlin), 
whereas two individuals (Mali Zdenci) were heterozygous 
with the a/c allele combination, and one with the a/c al¬ 
lele combination. At the Pgm-2 locus, 9 hybrid individu¬ 
als (Grubisno polje, Ulovcev mlin) expressed b/b alleles 
typical for P. lessonae, while two individuals (Grubisno 
polje, Ulovcev mlin) expressed the P. ridibundus-like vari¬ 
ant (c/c). A single introgression was recorded in the Aat 
locus at P. esculentus (Mali Zdenci) with the a/a pattern 
typical for P. ridibundus (Table 2). No gene introgression 
from P. ridibundus into P. lessonae and vice versa was found 
among the loci. 

Hemiclonal reproduction 

Pelophylax esculentus specimens were examined for allelic 
activity in testes and ovaries for three loci to estimate a ba¬ 
sic type of a gamete production in hybridogenetic hybrids. 
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Table 3. The expression of alleles for investigated loci in gonads of Pelophylax esculentus. Abbreviations: F - females, M - males, Aat - 
aspartate aminotransferase, Pgm-2 - phosphoglucomutase, Ldh-1 - lactate dehydrogenase; underlined alleles typical for P. ridibundus\ 
non-underlined alleles typical for P. lessonae. More intensely expressed alleles are in bold. 


Locality 

Ulovcev mlin 

Mali Zdenci 

Grubisno polje 

Allele 

F (n=7) 

M (n=3) 

F (n=l) M (n=6) 

F (n=9) 

M (n=14) 

Aat 

5 aa 

2 ab 

4 ab 

5 aa 

10 ab 


1 ab 



2 ab 


Pgm-2 

3 be 

2 bb 

2 be 

1 be 4 be 

4 be 

3 cc 

7 be 

3 bb 


1 cc 





Ldh-1 

7 cc 

2 cd 

1 ac 1 cd 

6 cc 

3 cd 




1 cd 

1 ad 

2 cd 




2 ad 


2 be 






2 cc 






2 ad 


Aat. Female hybrids showed a heterozygous somatic al¬ 
lele combination a/b except a single individual with the 
a/a allele combination from Mali Zdenci. Females from 
Grubisno polje and Ulovcev mlin translated the a/a allele 
combination in ovaries typical for P. ridibundus except for 
three of them that had the a/b allele combination. The in¬ 
dividual with the a/a allele somatic combination showed 
inactivity in oocytes for this locus. Males had somatic a/b 
allele combination with the same pattern in the testes. 

Pgm-2. Most hybrid females showed a heterozygous b/c 
pattern both in somatic tissues and ovaries, except for three 
females (Grubisno polje) expressing c/c profiles typical for 
P. ridibundus. The females with b/b or c/c allele combina¬ 
tion (Ulovcev mlin) in the somatic tissue displayed the 
same pattern in the ovaries. All hybrid males expressing 
b/c alleles in the somatic cells showed the same expression 
in their testes (Grubisno polje, Ulovcev mlin, Mali Zden¬ 
ci). Likewise, the males from Grubisno polje showing a so¬ 
matic b/b combination displayed the same allelic pattern 
in the gonads. 

Ldh-i. Females with somatic c/c allele combination ex¬ 
pressed the c/c pattern in gonads, as typical for P. ridibun¬ 
dus genome (Grubisno polje, Ulovcev mlin). A single indi¬ 
vidual shared a heterozygous a/c allele combination both 
in somatic tissues and gonads (Mali Zdenci). One hybrid 
female from Grubisno polje with the a/d allele somatic 
combination showed the same expression in ovaries. There 
was also one female from the same locality showing a b/c 
allele somatic expression, while the only c allele was visible 
in the gonads. 

The hybrid males with the c/d allele somatic combi¬ 
nation showed c/d (Grubisno polje, Mali Zdenci) or c/c 
(Grubisno polje) allele combinations in testes. Two males 
from Mali Zdenci with the somatic combination a/c and 
a/d showed the a/d pattern in gonads. Two males from 
Grubisno polje showing the b/c somatic allele combination 
expressed the b/c profile in their gonads. AH males express¬ 
ing c/d (Ulovcev mlin, Mali Zdenci, Grubisno polje) or a/d 
(Grubisno polje) somatic allele combinations had the same 
variant in testes. 


Most P. esculentus testes (89.8%) had both parental al¬ 
leles in their profiles, although P. ridibundus-specific allele 
was usually more intense; in total 6.1% had the P. lessonae- 
like pattern and 4.1% showed the P. ridibundus-like pattern. 
The most common pattern in the ovaries of hybrid females 
was that of the P. ridibundus-like type (66.6%). In 4.8%, the 
P. lessonae-like genome was expressed, while other gonads 
showed the pattern typical for hybrids (28.6%). Still, those 
alleles coming from P. ridibundus were more intense. Data 
are summarized in Table 3. 

Age structure and annual growth rate 

The age structure of the studied water frogs is summarized 
in Figure 2, whereas Figures 3 and 4 show a distribution 
of the estimated age and annual growth rate, respectively. 
The majority of P. ridibundus individuals (66.6%) were at 
least seven years old. In contrast, 90% of the P. lessonae and 
P. esculentus were placed into the age group between one 



Figure 2. The age structure of the investigated Pelophylax indi¬ 
viduals in Croatia. 
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Table 4. Summary of age estimates and number of lines of arrested growth (LAGs) of Pelophylax individuals. Abbreviation: n - number 
of individuals. 


Estimated age / Number of LAGs 

P. esculentus / n 

P. ridibundus / n 

P. lessonae / n 

2 

2 

0 

1 

3 

6 

1 

5 

4 

14 

0 

4 

5 

16 

0 

3 

6 

5 

0 

3 

7 

7 

1 

2 

8 

0 

1 

2 

9 

3 

1 

0 

10 

1 

0 

0 

13 

0 

1 

0 

Average age (mean ± standard error) 

5.1 ± 0.24 

8.0 ± 1.61 

4.8 ± 0.40 

Standard deviation 

1.75 

3.60 

1.79 

Total number 

54 

5 

20 



a Median 
□ 25 %- 75 % 
X Min-Max 



° Median 
□ 25 %- 75 % 
I Min-Max 


Figure 3. Box plots with the distribution of estimated age of the 
studied Pelophylax individuals in Croatia. 


Figure 4. Box plots with the annual growth rate (total body 
length/estimated age ratio) distribution for the studied Pelophy¬ 
lax populations in Croatia. 


and seven years. The youngest individuals were identified 
as being two years old {P. lessonae and P. esculentus), while 
the oldest frog was a P. ridibundus individual that had 13 
years. The oldest P. esculentus and P. lessonae individuals 
were estimated to have 10 years and eight years, respective¬ 
ly. The average age was 8.0 ± 1.61 years for P. ridibundus 
(n = 5; Table 4), 4.8 ± 0.40 years for P. lessonae (n = 20), and 
5.1 ± 0.24 years for P. esculentus (n = 54). We did not find a 
statistically significant difference between taxa in their age 
(n = 79, p = 0.096; Kruskal-Wallis test), nor in their annual 
growth rate (n = 79, p = 0.47). 

Discussion 

Pelophylax water frogs can easily move from one pond to 
another at least within 2.5 km distance (Juszczyk 1953). 


Distances between the studied localities are generally small 
and, therefore, it is likely that we studied a single popula¬ 
tion rather than several isolated populations. Genetic anal¬ 
yses supported the presence of all three taxa, P. lessonae, 
P. esculentus and P. ridibundus, so we classify the popula¬ 
tion as the R-E-L system. 

As in our case, hybrids usually dominate in most of the 
investigated R-E-L systems, likely due to their better vital¬ 
ity in the Fi hybrid state linked to a wider ecological tol¬ 
erance when compared to the parental species (Kotlik 
& SuLOVA 1994, Spasic-Boskovic et al. 1999, Borkin et 
al. 2002, Mayer et al. 2013). The representation of the in¬ 
vestigated population is described as follows: P. esculentus: 
P. lessonae: P. ridibundus - 9:3:1. 

The staining intensity of bands was asymmetrical in hy¬ 
brid specimens for the Ldh-i locus only. As we do not have 
other data measuring ploidy level, we cannot conclude yet 
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that the gene dose pattern observed is an effect of poly¬ 
ploidy. 

Although distinguishing between incomplete line¬ 
age sorting and genomic introgressions is a difficult task 
(Choleva et al. 2014), the episodic appearance of hybrid 
individuals with the recombinant genotype associated 
with observations of lessonae genes in the ridibundus ge¬ 
nome (UzzELL & Berger 1975, Uzzell et al. 1977, Tunner 
1979, Mezhzherin & Morozov-Leonov 1997) or rarely 
vice versa (Gunther & Hahnel 1976, Kotlik 1996), is at¬ 
tributed to the role of introgressions (Uzzell et al. 1977). 
Our observations of some introgressed loci in P. esculen¬ 
tus (Table 2), together with the presence of both parental 
genomes still present in the gonads of some hybrids (Ta¬ 
ble 3) suggest that hybridogenetic genome exclusion can 
be incomplete, take place in meiosis itself or be absent at 
all (Dolezalkova et al. 2016). As a result, the presence of 
both parental genomes in meiosis can be associated with a 
recombination followed by the production of at least partly 
recombinant gametes (Uzzell et al. 1977, Tunner & Hep- 
piCH 1981). Parallel explanation for the presence of both 
parental genomes in gonads at some loci may be a back¬ 
ground expression of somatic tissues forming analysed re¬ 
productive organs. At this point, crossing experiments or 
a detailed study of gamete production will need to be per¬ 
formed in order to understand the mechanism and type of 
gamete production in the Croatian population. 

Water frog populations from the L-E system normally 
contain both sexes of P. esculentus and P. lessonae. Popula¬ 
tions belonging to the R-E system mostly include hybrids 
of a male sex only (Graf & Polls-Pelaz 1989, Ragghian- 
Ti et. al. 2007). In natural populations throughout most of 
central Europe and northern Italy in general, P. esculentus 
females somewhat outnumber the males (Berger 1988). 
An explanation for this phenomenon could be the higher 
reproductive success of P. esculentus females. In fact, mat¬ 
ings between P. lessonae males and P. esculentus females re¬ 
sult in offspring with an even sex ratio. Conversely, crosses 
between P. lessonae females and P. esculentus males pro¬ 
duce an all-female progeny as a result of complex mech¬ 
anisms of sex determination, hybridogenesis and mating 
behaviour (Som & Reyer 2006). Contrary to our expec¬ 
tations to find more hybrid females, we found a larger 
number of males (23:17), which is a rather rare phenom¬ 
enon (Gunther et al. 1979, Tunner 1979, Wijnands 1979, 
Berger 1988, Spasic-Boskovic et. al 1999, Krizmanic & 
IvANOvic 2010). Our results also differ from populations 
found in the Pannonian Basin in Croatia, where a high 
number of females was observed by Berger (1988). How¬ 
ever, it should be noted that we collected our specimens 
during the breeding season, when water frogs usually show 
a strongly skewed operational sex ratio (OSR) with males 
outnumbering females (Wells 1977). Alternative explana¬ 
tions for our observation hypothesize that the genetic basis 
of sex determination in water frogs could be more complex 
than a simple XX-XY mechanism (Berger 1988). 

Despite only seven P. ridibundus were found, we con¬ 
sider the species as a self-perpetuating taxon because both 


males and females were found in the R-E-L population. 
The pattern of P. esculentus gamete production remains yet 
to be studied in detail; however a dominant allozyme pat¬ 
tern in gonads of hybrids was that of ridibundus-like. This 
suggests that P. esculentus produce ridibundus gametes and 
use lessonae gametes either from P. lessonae or from other 
hybrid for self-reproduction. Therefore, from a reproduc¬ 
tive perspective the R-E-L system is functionally a typi¬ 
cal L-E system. Although P. ridibundus can be produced 
by homotypic matings between two P. esculentus individu¬ 
als, the results would be female sex only (Hotz et al. 1992), 
which is not the case for this population. 

Data related to the age structure of water frogs are gen¬ 
erally rare. By comparing all three taxa within the R-E-L 
system our study results showed that the estimated aver¬ 
age age of P. ridibundus was relatively high (eight years), 
whereas P. lessonae and P. esculentus were between two and 
five years old. While the comparative data are not known 
for P. lessonae, P. ridibundus age estimates from Greece and 
Turkey ranged between 2.96 and 5.58 years (Yilmaz et al. 
2005, Kyriakopoulou-Sklavonounou et al. 2008, Gul et 
al. 2011), and P. ridibundus with P. esculentus from Poland 
were between 3.7 and 4.4 years of age (Socha & Ogielska 
2010). 

These data are close to our age estimates for P. lessonae 
and P. esculentus but not for P. ridibundus species (Fig. 2 
and 3). Currently, we have no clear explanation for the high 
age of studied P. ridibundus individuals. Despite that, our 
statistical analysis did not reveal differences in the estimat¬ 
ed age or in the annual growth rates between these three 
species. Also, we must note that our sample size of P. ridi¬ 
bundus was rather small and could have affected our statis¬ 
tical analyses. Thus, to get better idea about the relationship 
between the taxa and their age structure in the R-E-L sys¬ 
tem, more studies are needed that focus on species micro¬ 
habitat preferences, analysis of food competition and mat¬ 
ing success rates. 

The studied area represents one of the southernmost dis¬ 
tributions of the hybrid taxon P. esculentus not only in the 
Balkan Peninsula but considering its European distribution 
in general. The closest sites of the southernmost distribu¬ 
tion are known from Italy, Slovenia and Romania (Plot- 
NER et al. 2008). Data from Serbia (Krizmanic & Ivanovic 
2010) also help to define a southern line of P. esculentus' 
occurrence. Considering the active role of P. esculentus for 
interspecies transfer and the presence of a contact zone be¬ 
tween the southern limits of P. esculentus and the north¬ 
ern limits of P. ridibundus populations that are considered 
as still non-introgressed at least in mitochondrial DNA 
(Plotner et al. 2008), the region certainly deserves further 
attention in the context of studying P. ridibundus genetic 
resistance against P. esculentus-mediated introgressions. 
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